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Abstract. In this paper, we present mean gas and dust opacities relevant to the physical conditions typical of protoplanetary 
discs. As the principal absorber for temperatures below ~ 1 500 K, we consider spherical and aggregate dust particles of various 
sizes, chemical structure, and porosity, consisting of ice, organics, troilite, silicates, and iron. For higher temperatures, ions, 
atoms, molecules, and electrons are included as the main opacity sources. Rosseland and Planck mean opacities are calculated 
for temperatures between 5 K and 10 000 K and gas densities ranging from 10""* g cm"^^ to 10"' g cm"^ . The dependence on the 
adopted model of dust grains is investigated. We compare our results with recent opacity tables and show how different opacity 
models affect the calculated hydrodynamical structure of accretion discs. 
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1. Introduction 

Recently, significant progress toward the understanding of the 
possible composition and properties of dust grains and gas 
species in many astrophysical environments has been achieved. 
For instance, (sub)millimetre observations of molecular lines 
provided basic information about chemical composition and 
dynamical properties of the gas in discs around pre-main- 
sequence stars and young stellar objects (e.g., Bujarrabal et 
al. lT^ Olofsson, Liseau, & Brandeker l2IMI Thi et al. 2001 
Pietu et al. .2003,) . The infrared-to-millimetre continuum ob- 
servations of such environments constrain the properties of 
dust grains and can be used to estimate masses and thermal 
structure of the objects (e.g., Bouwman et al. .2000: Boogert, 



Hogerheijde, & Blake 2002 Tuthill et al.|2002|l. Finally, ex- 
perimental studies on the formation and spectra of various 



gaseous species (e.g., Butler et al. 120011 Sanz, McCarthy, & 
Thaddeus I2002> as well as the composition and properties 
of meteoritic, cometary, and interplanetary dust together with 
their laboratory analogues (e.g., Chihara et al. 2002 Mutschke 
et al. l2002> form a basis for theoretical investigations. 

On the other hand, the increase of computer power and 
new numerical methods stimulate the development of more 
sophisticated hydrodynamical models (e.g., Kley, D'Angelo, 



& Henning lMJTl D'Angelo, Henning, & Klev EUni Struck, 
Cohanim, & Willson I2002I I. Many of such simulations need 
an accurate treatment of the energy transport within the dusty 
medium (see e.g., Klahr, Henning, & Kley 119991 Niccolini, 
Woitke, & Lopez '2003 ), which requires a detailed description 
of the radiative properties of matter. Consequently, the adopted 
opacity model is an important issue. 

In this paper, we deal with physical conditions typical for 
protostellar nebulae and protoplanetary discs around low-mass 
young stellar objects. Virtually everywhere within the medium 
dust grains are the main opacity source, as they absorb radi- 
ation much more efficiently compared to the gas and because 
the temperature in these regions is low enough to prevent their 
destruction. However, for hotter domains (T ^ 1 500 K), where 
even the most stable dust materials cannot survive, it is neces- 
sary to take absorption and scattering due to molecular species 
into account. 

Recently, several extensive models describing the prop- 
erties and evolution of dust grains in protostellar cores 
and protoplanetary discs were proposed by Henning & 



Stognienko J1996t . Schmitt, Henning, & Mucha J1997t . and 
Gail (200]]Enn3- 
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Henning & Stognienko J 19961 ) studied the influence of par- 
ticle aggregation on the dust opacity in the early evolutionary 
phases of protoplanetary discs. They concluded that distribu- 
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tion of iron within the particles affects their optical properties 
in a great respect. Schmitt et al. (I1997> for the first time in- 
vestigated coUisional coagulation of dust grains in protostellar 
and protoplanetary accretion discs coupled with hydrodynam- 
ical evolution of these objects. They reported significant alter- 
ation of the thermal disc structure caused by the modification 
of the opacity due to dust growth. Gail ( I2001ll2002t considered 
annealing and combustion processes leading to the destruction 
of silicate and carbon dust grains consistently with the evolu- 
tion of a steady-state accretion disc. They found that the mod- 
ification of the dust composition in the inner regions due to 
these processes and its consequent transport toward outer disc 
domains aff'ect the opacity and, eventually, the entire disc struc- 
ture. 

A number of papers deal with the computation of 
Rosseland or/and Planck mean gas opacities in atmospheres 
of cool stars, protostars, and stellar winds. Alexander & 
Ferguson J1994t computed a set of opacity tables for temper- 
atures between 700 K and 12 500 K for several compositions. 
They considered a condensation sequence' of refractory mate- 
rials based on chemical equilibrium calculations and took into 
account absorption and scattering properties of these solids as 
well as various gas species. Note that the actual formation pro- 
cess is not modelled by such calculations. Helling et al. (2000 1 
calculated gas mean opacities for wide ranges of density, tem- 
perature and various chemical compositions based on up-to- 
date spectral line lists of the Copenhagen SCAN data base and 
studied the importance of the molecular opacity for the dynam- 
ics of the stellar winds of cool pulsating stars. An extension 
of this work will be used to construct our opacity model for 
protoplanetary accretion discs. 

In Appendix A, we give a brief overview of the most com- 
mon opacity models and studies and where they have been ap- 
plied. It can be clearly seen that there is a lack of papers which 
focus on calculations of both Rosseland and Planck mean opac- 
ities of grain and gas species for temperatures between several 
K and few thousands K in a wide range of densities based on 
both the best estimates of the dust composition and properties 
and recent improvements in molecular line lists. 

The first goal of this paper is to define such a model. In 
addition, we study how Rosseland and Planck mean opacities 
depend on the properties of dust grains and compare them with 
other opacity tables. Furthermore, we investigate how different 
opacity models affect the hydrodynamical structure of accre- 
tion discs. Our opacity model is freely available in the Internet: 
http://www.astro.uni-jena.de/Labomtory/labindex.html 

The paper is organised as follows. We introduce the opacity 
model in Sect. |2] The influence of the grain properties on re- 
sulting Rosseland mean opacities is described in Sect. 13. II The 
Rosseland and Planck mean opacities are compared to other re- 
cent opacity tables in Sect. 13. 21 We study how different opacity 
tables affect the hydrodynamical structure of active accretion 
discs in Sect. 13.31 A summary follows in Sect. |4] Finally, an 
overview of the opacity models cited in the paper is given in 
Appendix A. 

' which is not to be confused with a formation sequence (see also 
Woitke.2000j 



2. The opacity model 

2.1. Dust opacities 

In order to calculate frequency-dependent, Rosseland and 
Planck mean opacities of dust grains, we partly follow the 
schemes, proposed in Pollack et al. ( 119941 hereafter PHB), 
Henning & Stognienko ( II996I HS), and Voshchinnikov & 
Mathis 11999. VM). 

2.1 .1 . Composition model and grain sizes 

In this paper we follow the dust composition model for accre- 
tion discs by Pollack et al. (119941 Sect. 2.3) which is based on 
an analysis of a wide range of theoretical, laboratory, and obser- 
vational dust data. This composite model has been frequently 
used, for instance, in the evolutionary modelling of accretion 
discs around young stellar objects or estimates on their mass 
from millimetre observations (e.g.. Greaves et al. 1998, Jura & 
Werner [T999l D'Alessio, Calvet, & Hartmann l20UTl . 

The main dust constituents include amorphous pyrox- 
ene ([Fe, MgJSiOs), olivine ([Fe,Mg]2Si04), volatile and re- 
fractory organics (CHON material), amorphous water ice, 
troilite (FeS), and iron-. Following HS, we vary the rela- 
tive iron content in the silicates considering "iron-rich" (IRS) 
sihcates with Fe/(Fe-i-Mg)=0.4, "normal" silicates (NRM) 
with Fe/(FeH-Mg)=0.3, and "iron-poor" (IPS) silicates with 
Fe/(Fe+Mg)=0. However, the absolute amount of metallic iron 
in all these models is kept constant, which leads to the absence 
of solid iron in the first case and enhanced mass fraction of Fe 
in the third case. Such a variety of silicate models allows us to 
study the influence of iron content within the grain constituents 
on the extinction properties of dust. Another reason is that the 
exact mineralogical composition of the silicates in the proto- 
stellar clouds and protoplanetary discs is poorly constrained 
and can be different for various environments. Compared to 
HS, we re-estimated the absolute abundances of the silicates 
(FejMgi-j^SiOs, Fe2xMg2A-2Si04), iron, and troilite (FeS) in 
the cases of the "iron-poor" and "iron-rich" models from the 
Fe-stocheometric fractions keeping constant the total amount 
of iron. For the mass fractions of all dust constituents and their 
densities we follow Table 2 in PHB. Note the difference be- 
tween the iron mass fractions in the different dust models. 

The sublimation temperatures of the grain constituents are 
adopted from PHB (see Table 3 therein). We suppose that de- 
struction of dust materials occurs in a narrow range of temper- 
atures (~ 10 - 30 K). Given that the evaporation of the sili- 
cates and iron happens at approximately the same conditions, 
we do not distinguish between their evaporation temperatures 
and assume that they evaporate in one wide temperature range, 
~ 100 K. Thus, we account for six principal temperature re- 
gions: 

1. r ^ 155 K - all dust constituents are present; 

2. ~ 165K < r < 270 K - no ice; 

3. ~ 280K <T ^ 410 K - no ice and volatile organics; 
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4. ~ 440K < r ^ 675 K - silicates, iron, and troilite are 
present; 

5. ~ 685K < r $ 1 500 K - silicates and iron are present; 

6. ~ 1 500K < r < 10 000 K - gas-dominated opacities; 

We take into account the dependence of the evaporation 
temperatures of ice, silicates, and iron on gas density (the val- 
ues shown above are given for a density ~ 10"'" g cm"-'). Note 
that if one considers a condensation sequence of these materials 
(hysteresis behaviour), it results in higher values of the conden- 
sation temperatures compared to the evaporation temperatures. 

We assume that for the fifth temperature region the absolute 
amount of solid iron raises due to the destruction of troilite. The 
corresponding mass fractions of metallic iron are 6.15 x 10 ^*, 
2.42 X 10 and 1 .29 x 10"^ for the NRM, IRS, and IPS silicate 
models, respectively. 

As for the size distribution, we apply a modified MRN 
(Mathis et al. 1977) function suggested by Pollack, McKay, & 
Christofferson ( 1985 1. The modification consists of the inclu- 
sion of large (0.5//m - 5yum) dust grains. Such particle growth 
is expected to proceed efficiently at the early phases of the pro- 
toplanetary disc evolution due to the coagulation of small dust 
grains (see, for instance, Mizuno, Markiewicz, & VoeIk[l988 l. 
We do not consider other size distributions since the overall ef- 
fect of the particle sizes on the dust opacity is well studied (e.g.. 
Pollack et al. lT985l Beckwith, Henning, & Nakagawa^OOOl. 
It further allows us to compare directly our results with other 
works. 

2.1 .2. Grain structure and topology 

It becomes evident from theoretical investigations and labo- 
ratory experiments that the dust agglomeration is an efficient 
process in dense and relatively cold environments, like pro- 
tostellar cores or protoplanetary discs (e.g., Kesselman 119801 
Nuth & Berg 1994 Ossenkopf & Henning T^^T Kempf, 
Pfalz ner, & Henning 1999. Wurm & Blum 1998 2000, Blum 
et al. l2002> . Agglomeration leads to the formation of irregular 
particles consisting of hundreds or thousands of tiny subgrains. 
Usually, dust aggregates of two extreme kinds are considered, 
depending on the assumed coagulation processes, namely, PCA 
(particle-cluster aggregation) and CCA (cluster-cluster aggre- 
gation). As the laboratory and theoretical studies reveal, the 
PCA aggregates are sphere-like particles having a fractal di- 
mension of about 3. They have a compact "core" and a more 
rarefied "mantle". The CCA process results in the formation 
of very filamentary grains with complicated structure. They 
have fractal dimension of roughly 2 (Stognienko, Henning, & 
Ossenkopf 1995 1. 

During the evolution of parent objects, like molecular 
clouds or protostellar discs, chemical and physical processes 
can further modify the properties of dust grains. For in- 
stance, accretion of volatile materials on dust surfaces and 
their subsequent chemical processing are efficient in outer re- 
gions of protoplanetary discs and in protostellar clouds (e.g. 
Greenberg [T9671 Brown, Charnley, & Millar [T988l Hartquist & 
Williams 1990' Hasegawa & Herbst'1993", Willacy, Rawlings, 
& Williams .1994. Aikawa et al. ■1999j . This results in well- 



defined "core-mantle" or, more probably, "onion-like" grain 
structure. In protostellar discs, dust can be transported by the 
accretion flow toward hotter regions, where their volatile man- 
tle materials evaporate, and sputtering, annealing, combustion, 
and crystallisation processes may change the structure, com- 
position, and sizes of the grains (Gail 2001 2002i. Therefore, 
it seems obvious that the real astronomical grains must have a 
very complicated structure and topology. 

Unfortunately, modern computational methods and facili- 
ties allow only the consideration of somewhat simplified (but 
still reasonable) types of dust grains. In the present study, we 
focus on the following particle types: 

1. Homogeneous and composite aggregates; 

2. Homogeneous, composite, and porous composite spherical 
particles; 

3. Multishell and porous multishell spherical particles. 

An aggregate dust particle is assumed to be a cluster of 
small spherical subgrains sticking together. A particle is called 
"homogeneous" if it consists of only one dust component. On 
the contrary, "composite" means that a particle incorporates a 
fine mixture of various materials (heterogeneous particle). In 
addition to these two extreme cases of the chemical dust struc- 
ture, we consider a model of "multishell" grains, where each 
particle includes aU constituents distributed within concentric 
spherical shells. To study the influence of the porosity on the 
extinction properties of dust grains, we fill composite and mul- 
tishell spherical particles with vacuum. It is reasonable to as- 
sume that the optical behaviour of these porous multishell and 
porous composite particles may resemble that of more realistic 
kinds of dust grains. 

2.1.3. Computational methods 

The aggregate model and the numerical method to compute the 
optical properties of coagulated particles are adopted from HS. 
The aggregates are assumed to be in the form of PCA (50%) 
and CCA (50%) particles consisting of O.Ol/zm spherical sub- 
grains. The spectral representation of inhomogeneous media 
(Bergman ll978> and the discrete multipole method (DMM) are 
elaborated to calculate effective dielectric functions of aggre- 
gates (Stognienko et al. .l995.) . At first, we use the DMM to cal- 
culate the spectral function of the aggregated particles of a spe- 
cial topology, when their subgrains touch each other only at one 
point. Then, we account for the interaction strength between the 
subgrains (percolation), which varies with the size and type of 
the aggregates, by an analytical expression (see Equation 13 in 
HS). Finally, the effective dielectric functions of the aggregates 
are obtained by the spectral representation (HS, Eq. 5). In the 
case of the composite aggregates, we compute the optical con- 
stants of the composite material by the Bruggeman effective 
medium theory (Bruggeman ll935t . generalised to many com- 
ponents by Ossenkopf ( 195Tl. The optical properties of the dust 
aggregates are calculated with the usual Mie theory. It should 
be noted that this numerical approach is valid in the static limit 
only, which means that the scale of inhomogeneities within the 
particles must be small compared to the wavelength. Given the 
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Table 1. Mass fractions fj and densities pj of dust constituents 
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0.01/im size of the subgrains and the shortest wavelength of 
Q.lfim we considered, this condition is fulfilled. 

We use the numerical approach of VM to model (porous) 
composite and (porous) multishell dust particles. In this 
method, the composite grains are represented as spheres with 
many concentric shells, where each shell includes several lay- 
ers of randomly distributed dust materials. The multishell 
grains are modelled exactly as the composite ones but each 
shell includes only one layer of a dust constituent. Then a gen- 
eralised multilayered Mie theory can be applied to calculate 
their optical properties. As it has been shown by VM, a con- 
vergence in the optical behaviour of the multishell particles is 
achieved if the number of shells exceeds 3 and dust materi- 
als within shells are randomly ordered. In our calculations, we 
found that this number must be at least 20, because the highly 
absorbing material troilite is used, which induces interference 
within the shells and prevents fast convergence. Thus, in our 
case, a typical composite grain is modelled as a spherical par- 
ticle with about hundred shells. On the contrary, a multishell 
grain is represented by a spherical particle with only a few 
shells. 

We modify somewhat the dust model for the case of mul- 
tishell and composite spherical particles. We mixed the sil- 
icates and iron into one material with the Bruggeman mix- 
ing rule. A similar mixture of silicates, sulphides, and metals 
(GEMS, Glass with Embedded Metals and Sulphur) is found 
to be a common component of interplanetary dust particles 
(Rietmeijer & Nuth '2000). However, the main reason for this 
change is the convergence failure of the applied numerical 
method for the case of multishell grains with iron layers. This 
is due to numerical uncertainty which arises in the calcula- 
tion procedure for the Mie scattering coefficients in the case of 
highly absorbing materials, like iron (for further explanation, 
see Gurvich, Shiloah, & Kleiman 20011. 

We assume that each dust component has a total volume 
fraction in a particle according to its mass fraction and density, 
as specified in Table ^ For example, for the first temperature 
region and in the case of IPS silicate mineralogy, the mixture of 
iron and silicates occupies 8.9%, troilite - 1.6%, refractory or- 
ganics - 23.4%, volatile organics - 6%, and water ice - 60% of 
the entire particle volume, respectively. These values are sim- 
ilar in the case of NRM and IRS models. Thus, if the temper- 
ature is low, organics and ice are the dominant components of 
the dust grains. 



Unlike to the case of the composite particles, in the case of 
multishell spherical grains it is assumed that the distribution of 
dust materials is not random but follows their evaporation se- 
quence. Thus, for the first temperature region in the protoplane- 
tary disc, multishell spherical grains consist of a refractory core 
made of a mixture of silicates and iron and subsequent shells 
of troilite, refractory organics, volatile organics, and water ice. 
For higher temperature ranges, the number of shells is smaller 
since some materials are evaporated. In total, the number of 
shells in the case of multishell spherical particles varies from 
2 to 5. For the fifth temperature range (T > 700 K), where all 
troilite is converted to solid iron, we let this iron form an addi- 
tional layer on the grain surface. It is an extreme and probably 
physically unjustified case, but this allows us to study the in- 
fluence of the formation of a highly absorbing surface layer on 
resulting opacities. 

The porosity of particles is treated in a simple manner by 
the addition of 50% of vacuum (by volume) inside. For the case 
of the porous composite spheres, we consider vacuum to be one 
of the grain constituents, which forms additional "empty" lay- 
ers. In the case of the multishell grain model, we mix each shell 
with vacuum in the same way which is applied to create the 
composite spherical particles. That is, we subdivide each indi- 
vidual shell in many layers and fill some of them with vacuum 
according to the requested porosity degree. 

With the two computational approaches described above, 
we calculated the ensemble averaged absorption and scatter- 
ing cross section as well as albedo and the mean cosine of the 
scattering angle for all kinds of dust grains. Applying Eqs. 1- 
5 from Pollack et al. (1985 see Tabled of the present pa- 
per), the dust monochromatic opacity and, consequently, the 
Rosseland and Planck mean opacities were obtained for tem- 
peratures below roughly 1 500 K and a density range between 
10-'** gcm-^ and 10"^ gcm-^. A convenient analytical repre- 
sentation of the Rosseland and Planck mean opacities for every 
temperature region is provided as a 5 -order polynomial fit. This 
representation allows to calculate the Rosseland and Planck 
mean dust opacities accurately (~ 1%) and quickly for any 
given temperature and density values within the model applica- 
bility range, which is important for computationally extensive 
hydrodynamical simulations. The corresponding fit coefficients 
can be found in the code^. 
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2.2. Gas opacity 

The opacity of the very inner part of the protoplanetary accre- 
tion disc is dominated by various gaseous species. Here, the 
temperature is too high for dust to be present. 

Compared to the calculation of dust opacities, the calcula- 
tion of accurate Rosseland and Planck mean gas opacities is 
more challenging due to the large variation in frequency, tem- 
perature, and density of the absorption coefficient of numerous 
molecules, atoms, and ions. In addition, the body of data to be 
handled easily amounts to several millions of absorption lines 
per molecule. 

Missing data for absorption lines are critical for the cal- 
culation of Rosseland mean gas opacity since it is dominated 
by transparent spectral regions due to the harmonic nature of 
the averaging process. Therefore, each Rosseland mean is al- 
ways only a lower limit of the correct value. The opposite is 
true for the case of the Planck mean opacity - missing data for 
weak lines or bands cause an overestimation of the strong lines. 
Therefore, a Planck mean is always an upper limit of the case 
of ideally complete data. 

The dust opacity model for protoplanetary accretion discs 
outlined in the previous sections is supplied by a new ta- 
ble of gas opacities assembled on the basis of Helling ( 1999" 
Copenhagen SCAN data base) and Schnabel (2001, HITRAN 
data base). The gas opacity model is outlined in Helling 
et al. (I2000> and only a short summary is given here. The 
Rosseland and the Planck mean opacities are calculated from 
opacity sampled Unes lists. The data for the line absorption 
coefficients used in Helling et al. KOOO"! (CO - Goorvitch & 
Chackerian .1994, TiO - J0rgensen 1994; SiO - Langhoff & 
Bauschilder [T9931 H2 O - J0rgensen & Jensen [T9931 CH - 
J0rgensen et al. '1996; CN, C2 - J0rgensen & Larsson 119901 
C3 - J0rgensen 1989 HCN, C2H2 - J0rgensen 1990f were 
supplemented by data from the HITRAN 96 database (CH4, 
NH3, HNO3, H2CO, CO2, N2O, O3, SO2, NO2, HO2, H2, 
O2, NO, OH, N2). The opacity sampling of the latter was 
carried out by Schnabel J200II I. The set of continuum opac- 
ities and scattering includes continuum absorption from Hi 
(Karzas & Latte r 1961 1, (John 1988 1, H-hH (Doyle 1968 1, 
H ~ (Som erville fffpf . H+ (Mihalas ll965t . H e' (Ca rbon et 
al.[T569^, He i, C i, Mg i, Al i. Si i (all from PeachTWO^ as well 
as Thompson scattering on free electrons and Rayleigh scat- 
tering. Collision-induced absorption has been considered for 
H2-H2 and H2-He according to Borysow et al. (1997). 

The number densities of the ions, atoms and molecules are 
computed from an updated chemical equilibrium routine, in- 
cluding 14 elements and 155 molecules based on the JANAF 
table (electronic version of Chase et al. .1985 for more detail 
see Helling et al. I201j0l l. The element abundances are chosen 
mainly according to Anders & Grevesse ('1989'), but have been 
updated for various elements (see Helling et al. 2000). 

Using this approach, the Rosseland and Planck mean 
gas opacities were computed for temperature ranges be- 
tween 500 K and 10000 K and for gas densities between 
~ 10''^ gem'-' and ~ 10'^ gem'-*. In contrast to the dust 
opacity, no simplified analytical expression can be found for 
the gas mean opacities because of their sensitive dependence 



on temperature and density. Thus, we apply a second-order in- 
terpolation in order to calculate the gas opacities for any given 
temperature and density value from tabulated values. 

2.3. Opacity table 

In order to assemble the opacity table, we take into account ei- 
ther only dust opacity data for low temperatures, T ^ I 500 K, 
or only gas opacity data for higher temperatures. As has been 
shown by many authors (e.g., Lenzuni, Gail, & Henning ll995t . 
it is an accurate approach because dust dominates the ab- 
sorption properties of matter whenever it is present. It has 
been shown that already 30% of total SiO content in the solid 
phase is enough to exceed even the Planck mean gas opacity 
(Helling[l999j. This fraction will decrease with increasing im- 
purity of the grains. However, for the dust-to-gas transitional 
region (T from ~ 1 400 K to 1 600 K), where the last dust grain 
population evaporates, it is necessary to calculate the opacity of 
dust and gas simultaneously. In this narrow temperature range 
the resulting Rosseland and Planck mean opacities are certainly 
going down by few orders of magnitude, so one can apply a 
simple linear interpolation to estimate their values with a good 
accuracy. Still, we note that our table may give only approxi- 
mate opacity data for that temperature range. 



3. Results 

In this section, we present the opacities of all types of dust 
grains and discuss the differences between them. Second, we 
compare the Rosseland and Planck mean opacities with other 
recently published opacity models. Finally, we study the influ- 
ence of the adopted opacity model on the accretion disc struc- 
ture. 

3.1. Opacity and the dust models 

The Rosseland mean opacities /cr computed for all dust mod- 
els are presented in Fig.[T] We compared two silicate models, 
namely, the IPS (left panel) and the IRS (right panel). Shown 
are the following dust models: multishell spheres, composite 
aggregates, composite spheres, homogeneous aggregates, ho- 
mogeneous spheres, porous multishell grains, and porous com- 
posite spheres. 

The most prominent trends in behaviour of the Rosseland 
mean dust opacities can be summarised as follows: 

1 . There is a significant difference in the calculated dust opac- 
ity values between the aggregates, (porous) composite and 
(porous) multishell spherical particles for the first (T ^ 
150 K) and fifth (T <: 700 K) temperature regions; 

2. For intermediate temperatures, 150 K $ T $ 700 K, the 
resulting Rosseland mean opacities do not show profound 
variations with the applied dust models, except for the case 
of the IPS homogeneous aggregates; 

3. The discrepancy between opacity curves are smaller for the 
case of the IRS silicate model compared to the IPS; 
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Fig. 1. The Rosseland mean opacities in the case of IPS (left panel) and IRS (right panel) silicate dust model. The following 
dust particle types are indicated: multishell spheres - solid line, composite aggregates - triangles, composite spheres - circles, 
homogeneous aggregates - dashed line, homogeneous spheres - dot-dashed line, porous multishell grains - dotted line, and 
porous composite spheres - pluses. 



4. Addition of vacuum inside compact composite and multi- 
shell spherical particles leads to significant modification of 
their opacities. 

In the first temperature region, organics and water ice are 
the dominant dust materials according to our compositional 
model (see Table 0. At that temperatures {T ^ 150 K), the 
main contribution to the Rosseland mean opacity comes from 
monochromatic opacities in the wavelength range about 30/vm 
- 400//m. In this spectral range, iron and troiUte have higher 
values of the real (n^) and imaginary (A:^) indices of refraction 
compared to the other materials considered in our model. For 
instance, ai A - 100//m, iron has n^^j^^ = 95.02, k^^^Q^^ - 
181.95, troilite has n\^^^^^ = 8.5, k\^^^ = 0.73, whereas or- 
ganics and water ice have k^qq^^ - 2.14, ^fo'o/im ~ ^-^^ '^^'^ 
"loVm ~ ^-^2, ^{oo^j^ - 0.05, respectively. Furthermore, for 
troilite and especially for iron these values rapidly increase 
with wavelength. Thus, the resulting optical properties of a dust 
grain become sensitive to the absolute amount of Fe and FeS 
and its distribution inside the particle. In such a case, different 
theories to calculate optical properties of an aggregate particle 
may give different results and should be used with caution (see 
Stognienko et al. [T995l Fig. 6, 8). 

We considered aggregated particles of two kinds, namely, 
composite and homogeneous aggregates. The optical constants 
of the composite material do not show a peculiar behaviour at 
long wavelengths. Moreover, the refractive index of this com- 
posite changes only little if one is switching from the IPS to the 
IRS siUcate model. One reason is that the amount of solid Fe is 
small compared to other constituent materials. Another reason 
is that the total iron abundance is kept constant in all compo- 
sitional models. Hence the resulting optical properties of the 
composite aggregates are not very sensitive to the actual topol- 
ogy of the particles and adopted silicate model. Indeed, as it 
can be clearly seen in the Figure, the dust opacity values of the 
composite aggregates (triangles) are close to the a-r of the ho- 
mogeneous spheres (dot-dashed line). The maximum deviation 



of these opacity curves is achieved in the case of the IPS sili- 
cates for temperatures higher than ~ 700 K. Here, the absolute 
amount of metallic iron is increased due to conversion of FeS 
to Fe. In addition, the Rosseland mean opacities of the com- 
posite aggregates do not depend much on the adopted silicate 
model (compare triangles on the left and right panel). 

On the contrary, the a-r values for the homogeneous ag- 
gregates do demonstrate a strong variations with the compo- 
sitional model. In the case of the "iron-rich" silicates (right 
panel) the Rosseland mean values of the homogeneous aggre- 
gates (dashed line) he much closer to the kr of the composite 
aggregates (triangles) then for the "iron-poor" compositional 
model (left panel). As it has been shown already by HS, this is 
caused by the presence of bare iron aggregates in the case of the 
IPS homogeneous aggregate model. Due to the extremely high 
absorptivity induced by the strong interactions between indi- 
vidual aggregate sub-grains, the optical properties of such iron 
clusters determine the overall behaviour of the resulting opac- 
ities. Note, that in the case of the "iron-rich" silicate composi- 
tion, all iron is locked inside silicates and the absolute amount 
of troilite is reduced by a factor of two compared to the IPS 
and NRM models. In the absence of a population of highly 
absorbing grains, the Rosseland mean values of the homoge- 
neous aggregates are close to that of the composite aggregates. 
However, at temperatures higher than ~ 700 K, troilite is con- 
verted to solid iron which increases the ar values of the homo- 
geneous aggregates in respect to that of the composites (com- 
pare dashed lines with triangles on the right panel of Fig.QJ. 

The situation is different for the case of the (porous) com- 
posite and (porous) multishell spherical particles. As it was 
mentioned in Sect. 12.1.31 we changed the compositional model 
in this case and locked all solid iron in silicates using the 
Bruggeman rule of the effective medium theory. Thus, the only 
material with a high absorptivity at long wavelengths which re- 
mains in all compositional models is troilite. The metallic iron 
is another highly-absorbing dust component, but it is present in 
the fifth temperature region only, at J > 700 K. 
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Voshchinnikov & Mathis J1999> have shown that conven- 
tional EMTs are rather approximate in the case of small com- 
posite spherical particles, x - InajA «; 1, when one of the 
grain constituents has a large refractive index (see Fig. 3b 
therein). Given a typical size a - O.lyum of our dust grains, 
a typical wavelength of A - lOOyum for the first temperature re- 
gion, and the high refractive index of troilite at that wavelength, 
this condition is fulfilled. Thus, we adopted the approach of 
VM to model a composite grain as a spherical particle with 
many concentric layers and applied it for all temperature re- 
gions. Note that there is always an interference between the 
layers, which makes the optical properties of such composite 
particles to be different from that of composite grains with well 
mixed dust constituents. This is especially true if one of the 
dust components has a high absorptivity. In our compositional 
model for layered spherical particles, this is troilite in the first 
temperature region while in the fifth region it is iron. Therefore, 
we may expect to see the difference between the Rosseland 
mean opacity values of the composite spheres and composite 
aggregates, particularly for T ^ 150 K and T > 700 K. 

As it can be seen in Fig.Q] indeed the /cr values in the case 
of the composite spherical particles are much higher than for 
the composite aggregates. For instance, at T = 10 K this differ- 
ence can reach factors of 10 and 20 for the IPS and IRS com- 
positions, respectively (compare circles and triangles on the 
left and right panels). For higher temperatures, the Rosseland 
mean opacity curves of the composite spheres and aggregates 
are close to each other till T ~ 700 K is reached. In this temper- 
ature region, we assumed that iron forms a layer on the surface 
of the composite spherical grains. Such a layer "screens out" 
all underlying materials and totally dominates the optical be- 
haviour of the entire particle. Due to this fact, the dust opacity 



values of the composite spheres in the fifth temperature region 
is nearly the same for both the IPS and IRS compositions. 

The multishell spherical particles have a restricted number 
of layers compared to the composite spheres, namely from 2 
to 5, depending on the temperature region. The troilite layer 
is assumed to be the first layer after the silicate core and thus 
troilite is "hidden" inside. It prevents a strong interference be- 
tween the consequent particle shells as it is the case for the 
composite spheres. Then one may expect that the kr values of 
the multishell spheres should be lower than that of the compos- 
ite spherical particles, especially for the IPS model. As it can 
be clearly seen in Fig.^ this is true for T $ 150 K, whereas 
for higher temperatures both opacity curves almost coincide 
(compare solid line with circles). Hence, the actual distribution 
of dust constituents within a multilayered spherical particle is 
not that important for the relevant Rosseland mean opacities at 
T < 150 K. 

The addition of vacuum inside the compact composite and 
multishell spherical grains leads to a significant increase of 
the corresponding a-r values for the first and fifth tempera- 
ture regions (compare solid line with dotted line and circles 
with pluses). The first reason is that the density of the porous 
grains becomes lower than the density of the compact parti- 
cles. Second, for the porous spheres the relevant extinction ef- 
ficiencies are higher compared to that of the compact spherical 
particles if some of the dust constituents have a particularly 
high absorptivity, like troilite in the first and iron in the fifth 
temperature region, respectively. This is due to the coherence 
between the particle layers. The interference is more intense 
for the case of the composite sphere since it has more con- 
centric layers and a nearly homogeneous distribution of the 
dust constituents from the centre to the surface compared to 
the multishell spherical particle. Note that in the fifth tempera- 
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ture region, both composite and multishell spheres have a sim- 
ilar topology, namely, a silicate core covered by iron mantle. 
Therefore, it is naturally to expect that they have a similar be- 
haviour of the resulting Rosseland mean opacities. 

As it can be clearly seen in Fig. ^ the kr values of the 
porous composite spheres are higher than that of the compact 
composite spherical particles at T < 150 K and T > 700 K by 
a factor of 2 for both the IPS and IRS compositions (compare 
pluses and circles). For other temperatures, the corresponding 
Rosseland mean curves lie close to each other. The situation 
is similar for the case of the multishell and porous multishell 
spherical particles (compare solid line with dotted line). As we 
expected, in the fifth temperature region the Rosseland mean 
opacities for the case of the porous composite and porous mul- 
tishell spheres have almost the same values. 

3.2. Comparison to other studies 

In Fig. 12 the Rosseland (left panel) and Planck (right panel) 
mean opacities composed of the NRM composite aggregate 
for the low temperature range and gas opacities for the high 
temperature range are compared with other models. We plotted 
these values for a wide temperature range, T from ~ 10 K to 
10^ K and for gas densities which scale as 10"'^ x gcm"^. It 
allows us to make a comparison between the models in a wide 
temperature as well as density range simultaneously. 

As it is clearly seen in Fig. |2l (left panel), the discrepancy 
between the Rosseland mean values provided by various mod- 
els is negligible at high temperatures (T <; 3 000 K). On the 
contrary, the Planck mean opacity values differ by a few order 
of magnitude in this temperature range (compare curves on the 
right panel). The reason is that the Rosseland mean is much 
less sensitive to differences in the material data than the Planck 
mean due to the nature of the averaging process. However, the 
Planck mean heavily depends on the adopted values of the band 
and line strengths which vary for different line lists and on the 
adopted chemical equilibrium constants^ 

Therefore, the Kp of the Opacity Project (OP, crosses) 
are much larger than all the other opacity models since it 
combines more atomic opacity sources (see e.g. Table 3.3 in 
Helling 1999 1. This model does not contain molecules for 
temperatures 1 500 K - 5 000 K, which become important 
absorbers in this temperature range. The difference between 
the Planck mean opacity values in the case of the AF model 
(dashed line) and our model (solid line) may be caused by 
different molecular line data and a different set of chemi- 
cal species adopted. The same is true for the Kp values of 

Note that the chemical equilibrium constants, Kp, used by differ- 
ent authors can cause differences in the resuiting opacity values be- 
cause it affects the number density of species. The same effect will 
be caused by the neglect of the metal ions in the chemical equilib- 
rium calculations (see, e.g., Helling et al. '2000 1 . It is, however, not 
very straightforward to decide which line list is the most correct (for a 
discussion see J0rgensen 2003 i. We nevertheless dare to demonstrate 
the difficulties arising from comparing Planck mean opacities calcu- 
lated by different authors because it has ~ to our knowledge - not been 
pointed out clearly in the literature. Only comparisons of Rosseland 
mean opacities are presented, e.g. by Alexander & Ferguson cl994< . 



Sharp J 19921 diamonds), which lay somewhat in between the 
Planck mean opacities of our model (solid line) and the values 
of Alexander & Ferguson (dashed line). 

For temperatures lower than about 1 500 K, dust grains are 
the main opacity source. As it has been shown by Pollack et 
al. (1994 1, in this case the difference between the Rosseland 
and Planck mean opacities computed for the same model is 
small, ~ 30% (see Fig. 4b therein). The reason is that both 
opacities are dominated by continuum absorption and scatter- 
ing rather than absorption lines in this case. In what follows, we 
focus on the low-temperature Rosseland mean opacities only 
(T ^ 3 000 K). 

The Rosseland mean opacity values kr computed by the 
model of Bell & Lin (|][994 BL) strongly deviate from the 
A-R calculated by other models. For example, this difference 
can reach a factor of hundred for T ~ I 500 K - 1 800 K 
(compare dot-dashed and solid lines in Fig. The reason is 
that the BL data are a modification of the old analytical Lin 
& Papaloizou (.1980j opacities, which are based on the opac- 
ity tables of Alexander J1975> and Cox & Stewart J1970al 
I970bi supplemented by the data from Alexander, Augason, 
& Johnson f 1989 1 for T < 3 000 K. The model of Bell & Lin 
includes dust grains of two types, namely, homogeneous icy 
and metallic spherical grains. Since we used the more advanced 
dust compositional model of PHB and a different set of evapo- 
ration temperatures, this causes lower Kp, values in the case of 
the BL model compared to the other models for T ^ I 500 K. 
For higher temperatures, ~ 1 500 K - 3 000 K, this deviation 
is huge, ~ 100 times. As it has been shown by AF, the reason 
is that Bell & Lin truncated monochromatic opacities of water, 
which is one of the main absorbers at such temperatures, at a 
too short wavelength in their calculations. The missing opacity 
data affect the resulting Rosseland mean opacity. They do, so 
far, not draw any conclusions regarding Planck mean opacities. 

In overall, the opacity curves of all other models do not 
show such a strong difference between each other (compare 
triangles, dashed line, solid line, and line with open squares). 
Our dust model differs from the model of Pollack et al. by tak- 
ing into account an aggregate nature of cosmic dust grains and a 
new set of optical constants, but dust size distribution, composi- 
tion, and evaporation temperatures are the same. However, the 
difference between the Rosseland mean opacity values of these 
two model can reach about a factor of two (see also Fig. 5a in 
Henning & Stognienko ll996> . 

The /cr values of the opacity model of Alexander ("1975' line 
with open squares) is lower for the dust-dominated temperature 
region (T ^ 1 500 K) and higher for the dust-to-gas transi- 
tional region (1 500 K ^ T 5 1 700 K) compared to our model 
(solid line) by factors of 5 and 100, respectively. The reason 
is, as mentioned above, that we used an approximation to com- 
pute Kp in that temperature region, where the last dust grains 
get evaporated, which is not a very accurate approach. On the 
contrary, the opacity model of Alexander assumes the presence 
(but no evaporation!) of small 0.1/im spherical silicate grains 
in a rather approximate way, assuming that all dust is homo- 
geneously condensed when the gas becomes saturated. Since 
here a phase-transition takes place, a supersaturated gas would 
be needed which results in higher molecular abundances than 
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AF derive from their equilibrium consideration (for a discus- 
sion see e.g. Woitke & Helling 2003 1. Moreover, the conden- 
sation begins at higher temperature than the value we assumed 
for the evaporation of the last grain constituents because of the 
dust hysteresis. Finally, the model of AF neglects the presence 
of other refractory materials, like iron, in the dust-dominated 
temperature range, which makes the relevant Rosseland mean 
opacity values lower than provided by our model. 

On the contrary, the model of AF does consider several 
dust materials, namely, iron, silicate, carbon, SiC assumed to 
be present as small ellipsoidal grain particles. Therefore, the 
corresponding Rosseland mean dust opacities nearly coincide 
with our values (compare dashed and solid lines, respectively). 
The same is true even for the dust-to-gas transitional region, 
where our opacity model gives rather approximate opacity val- 
ues. 

3.3. Opacity and disc structure 

We compare the thermodynamical structure of a typical proto- 
planetary disc around a low-mass star computed with two dif- 
ferent opacity tables in Fig.|3] The 1 + ID model of an active 
steady-state accretion disc of Ilgner (2003 1 was used with the 
following input parameters: M* = IMq, M - 10'''Mq yr"', 
and a - 0.01. Here, M* is the stellar mass, M is the mass ac- 
cretion rate, and a is the parameter describing the kinematic 
viscosity. Note that this model incorporates a star as a gravita- 
tional center only and does not take into account the effect of 
the stellar radiation on the disc structure. 

The thermal structure shown on the left panel was obtained 
with the Rosseland mean opacity table of Bell & Lin (1994'). 
On the right panel, we present the same disc structure but for 
the case of the IPS homogeneous aggregate model (IPSHA) 
supplemented by the gas mean opacity. We choose these two 
opacity models as the overall difference between them is the 



largest among the different models (compare solid and dash- 
dotted lines in Fig.|5}- 

It can be clearly seen that the higher values of the Rosseland 
mean opacity in the case of the IPSHA model leads to a hot- 
ter and more extended disc structure. For instance, the vertical 
scale height of the disc at 20 AU is equal to 1.3 AU for the 
former and 2.1 AU for the latter opacity models, respectively. 
Consequently, there is also a variation of the density structure 
between the models, namely, the disc density is higher for the 
BL opacity model compared to that of the IPSHA model. The 
temperature difference is also prominent. For example, the mid- 
plane temperature of 100 K, which roughly corresponds to the 
ice melting point, is reached at 6.5 AU for the model of Bell & 
Lin, whereas in the case of the IPSHA model it is at w 8 AU. 

To confirm our findings, we did a similar comparison with 
another code. We used a full 2D hydrodynamical code designed 
to simulate the interaction of the protoplanetary disc with a pro- 
toplanet (D' Angelo 2001 ). The parameters of the model were 
as follows: = IM©, Md = O.OIM©, v = 10'^ cm^s^^, 
ju - 2.39, and -y - 1.4, where Md is the total disc mass, v is 
the kinematic viscosity, fj. is the mean atomic weight of the gas, 
and 7 is the adiabatic exponent. 

The Bell & Lin opacity model together with IPS and IRS 
homogeneous aggregate models were chosen for the compar- 
ison. The midplane temperature for all three opacity models 
is shown in Fig. |3 It can be clearly seen that the difference 
in the temperature values between all models can reach about 
50% for the disc radii rd ^ 2 AU. Note that it disappears at 
larger distances, ~ 10 AU. As expected, in the case of the 
BL model the midplane temperature is the lowest almost ev- 
erywhere (dashed line), whereas for the IPS homogeneous ag- 
gregates it is the highest (solid line) and the IRS temperature 
values lie between them (dotted line). The reason is the same 
as for the case of the 1 -H ID disc model, namely, lower opac- 
ity values of the BL model compared to both the IPS and IRS 
opacities and higher values of IPS opacity in comparison with 
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Fig. 4. The midplane temperature of the accretion disc obtained 
with the BL opacity model (dashed Hne) and IPS (soHd Hne) 
and IRS (dotted line) homogeneous aggregate dust models. 

those of BL and IRS models. It is interesting that at T ~ 130 K 
{r ~ 2 AU) all the temperature curves are very close to each 
other. This is due to the adopted ice evaporation temperatures, 
which is a little lower in our case compared to the model of 
Bell & Lin. That leads to nearly the same opacity values for all 
three models at a restricted set of temperatures around ~ 130 K 
(compare solid and dot-dashed lines in Fig.|3}- For r ^ 10 AU 
the temperature curves are close to each other because the cor- 
responding opacity models have rather similar low-temperature 
opacity values. 

Thus, we showed that the difference in the Rosseland opac- 
ity tables applied in hydrodynamical calculations leads to devi- 
ations in the resulting disc structure. As it has been shown by 
Markwick et al. (,2002.) . the temperature distribution is a crucial 
factor for the chemical evolution in the inner parts of accretion 
discs. Therefore, proper opacity modelling is an important is- 
sue in order to follow the dynamical and especially chemical 
evolution of protoplanetary accretion discs. 

4. Summary 

We have compiled the Rosseland and Planck mean dust and 
gas opacities for the temperature range T e [5, lO**] K and 
for gas densities p 6 [10""^, 10"^] gem"-' which is appropriate 
for the conditions in protoplanetary discs. The absorption and 
scattering due to dust grains of different compositions, shapes 
and topological structures and the absorption provided by over 
30 atoms, molecules, ions were taken into consideration. The 
corresponding well-documented numerical code together with 
representative data and figures are electronically available: 
.■//www.astro. uni-jena.de/Laboratory/labindex.html 
It has been found that the topological distribution of highly 
absorbing materials, such as iron and troilite, within dust 
grains dominates the relevant optical properties to a large de- 
gree. Particularly, the composite and porous composite spher- 
ical grains modelled in a special way have remarkably high 
Rosseland mean opacities at T 5 150 K and T ^ 700 K even in 
comparison with the composite aggregated particles. We have 
shown that at intermediate temperatures the Rosseland mean 
opacities of distinct grain models are close to each other It has 
been demonstrated that the difference between the opacity val- 



ues of various dust models is smaller in the case of the compo- 
sitional model with a smaller amount of solid iron and troilite. 
We found that porous composite and porous multishell spher- 
ical particles show higher opacity values in comparison with 
their compact analogues. 

We performed a comprehensive comparison of our results 
with other recent opacity models. We found a significant dif- 
ference between the opacity models in the case of the Planck 
mean and a good agreement between them for the case of the 
Rosseland mean at T ^ 1 500 K, where gas species are the 
main opacity sources. For lower temperatures, where opacities 
are dominated by dust grains, there is a discrepancy (a factor of 
~ 3 at most) in both the Rosseland and Planck mean values for 
all considered models. 

We demonstrated that differences in the Rosseland mean 
opacity values provided by distinct opacity models affect the 
hydrodynamical structure of active steady-state accretion discs. 
Namely, higher values of the Rosseland mean opacity lead to a 
hotter and more extended disc structure in the case of 1 + ID 
and 2D disc modelling. 
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Appendix A: An overview of tlie opacity models 

The goal of this compilation is to provide essential information 
about opacity models cited in the text. We show for which ele- 
mental composition they are developed, what kind of opacities 
they supply, and in what temperature and density ranges they 
work. In addition, key references to the studies, where these 
opacity models have been applied, are given and the aims of 
the investigations are briefly mentioned. 

Primarily, we distinguish between two kinds of opacity 
models. The first kind is designed for stellar evolutions, where 
it is more convenient to use a special parameter, R - p/T^,T(, - 
r/10* K instead of gas density p (see discussion in Rogers & 
Iglesias ,1992. 1. Thus, opacity data of such models are assem- 
bled in rectangular tables in R-T(, space. However, in hydro- 
dynamical simulations of accretion discs, it is more convenient 
to have opacity tables composed in p-T parameter space, for 
which we refer as to the second type of the models. The direct 
conversion of R to p in the opacity models of the first kind leads 
to trapezoidal tables, where opacity values for different temper- 
atures have different density intervals. We mark such models 
with "a" and show the maximum limits of density and temper- 
ature values for them. 

In addition, we point out if opacity data are available on- 
line in the Internet or via E-mail ("b"). By default, the pa- 
pers are supposed to contain the opacity data in a tabular form. 
Otherwise, when only opacity plots are available, we label the 
corresponding papers with "c". The analytical opacity models 
are marked with "d". 



mip 



D. Semenov et al.: Opacities for protoplanetary discs 



Table A.l. Opacity models cited in the text 
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" molecular and atomic opacities, 

" a trapezoidal table in p~T parameter space, 

it can be retrieved via E-mail or the Internet, 

only plots are available, 

analytical expressions are used, 
° LMC, SMC mean the element abundances of the Large and the Small Magellanic Clouds 
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